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MHD Sausage Mode Oscillations in Coronal Loops 
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ABSTRACT 

co 
o . 

A recent study by Nakariakov et al. pointed out that the dispersion relation of MHD sausage 
mode oscillations has been incorrectly applied to coronal loops, neglecting the highly dispersive 
0_i nature of the phase speed and the long- wavelength cutoff of the wave number. In the light 

of these new insights we revisit previous observations that have been interpreted in terms of 
MHD sausage mode oscillations in coronal loops and come to the following conclusions: (1) 
Fast sausage MHD mode oscillations require such a high electron density imposed by the wave 
number cutoff that they can only occur in flare loops; (2) In the previously reported radio 
observations [y w 100 MHz to 1 GHz) with periods of P « 0.5 — 5 s, the fast sausage MHD 
mode oscillation is likely to be confined to a small segment (corresponding to a high harmonic 
node) near the apex of the loop, rather than involving a global oscillation over the entire loop 
length. The recent microwave and soft X-ray observations of fast periods (P ~ 6 — 17 s) by Asai 
et al. and Melnikov et al., however, are consistent with fast sausage MHD oscillations at the 
fundamental harmonic. 
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% ■ 1. INTRODUCTION 

Coronal seismology became an efficient new tool that uses standing MHD waves and oscillations as 
a tool to explore the physical parameters of the solar corona. There are three basic branches of solutions 
of the dispersion relation for propagating and standing MHD waves: the slow mode branch (with acoustic 
phase speeds), the fast mode branch and the Alfven branch (with Alfvenic phase speeds). Furthermore, 
each branch has a symmetric and asymmetric solution, termed the sausage and kink mode (Roberts et 
al. 1984). All of these MHD oscillation modes have been now detected with imaging observations in 
recent years: transverse fast kink-mode oscillations with TRACE (Aschwanden et al. 1999; Nakariakov et 
al. 1999); longitudinal (slow-magnetoacoustic) modes with SUMER (Wang et al. 2002; Ofman & Wang 
2002), and fast sausage-mode oscillations probably with the Nobeyama Radioheliograph (Asai et al. 2001; 
Melnikov, Reznikova, & Shibasaki 2002). The latter type of fast sausage- mode oscillations is the least 
established one, although numerous non-imaging radio observations exist that have been interpreted earlier 
in terms of this mode. However, a recent study by Nakariakov et al. (2003) pointed out that the dispersion 
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relation and oscillation period has been incorrectly applied to the data, because the highly dispersive nature 
of the phase speed and the long-wavelength cutoff in the wave number has been ignored. Here we revisit 
the relevant earlier radio observations in the light of these new insights and attempt to derive the correct 
physical parameters. 



2. THEORY 



The dispersion relation for magnetosonic waves in cylindrical magnetic flux tubes has infinitely many 
types of long-wavelength solutions in the fast-mode branch (n = 0, 1, 2, 3, ...), with the lowest ones called 
sausage modes (n = 0) and kink modes (n = 1). The solutions are depicted in Fig. 1, showing that the kink 
mode solution extends all the way to the long-wavelength limit ka i— > 0, while the sausage mode has a cutoff 
at a phase speed of v p h =VAei which has no solution for small wave numbers ka < k c a. The propagation 
cutoff for the sausage mode occurs at the cutoff wave number k c (Edwin & Roberts 1983; Roberts et al. 
1984), 
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where j , s = (2.40,5.52, ...) are the zeros of the Bessel function J . Here we restrict our attention to the 
genuine sausage mode with transverse mode number s = 1, i.e. jo,i- The cutoff frequency u at the cutoff is 
k c VAe = u c . Under coronal conditions, the sound speed amounts typically to c s = 150 x \JTmk ~ 150 — 260 



km s 1 (with T« 1-3 MK), which is much smaller than the typical Alfven speed (va 
we have c s <C va- In this case the tube speed, 



2„,2 



C T = 



CZV 



1000 km s- 1 ) so 



(2) 



is similar to the sound speed, c T < c s , and the expression for the cutoff wave number k c (Eq. 1) simplifies to, 
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In the low-/? corona, the thermal pressure is much smaller than the magnetic pressure, so that we can 
assume almost identical magnetic field strengths inside and outside of coronal loops, i.e., B ~ B e , so that 
the ratio of external to the internal Alfven speed vac/va is essentially a density ratio ^Jn a /n e , 
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From this expression we see that for typical density ratios inferred in the solar corona, e.g., 
n e /no ~ 0.1 — 0.5 (Aschwanden et al. 2003), the cutoff wave number k c a falls into the range of 
0.8 < k c a < 2.4. Therefore, we would expect that a long-wavelength sausage mode oscillation is completely 
suppressed for the slender loops for which kink mode oscillations have been observed, which have wave 
numbers of ka = 2ira/\ = (n/2)(w/l) w 0.04— 0.08 (Aschwanden et al. 2002). Coronal conditions with 
h« 1 were also inferred from fast kink mode oscillations (Nakariakov & Ofman 2001). Of course, higher 
harmonics would have correspondingly shorter wavelengths, i.e. k = Nw/l with = 1,2,3, ...). On the 
other side, the short-wavelength limit is given by the maximum possible loop width (a = w/2 < I), so the 
maximum wave number is ka < k m a = 27r(a/A) m = n(a/l) m — ir. The occurrence of global sausage mode 
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oscillations therefore requires special conditions, k c < k < k m , i.e., very high density contrast no/n e and 
relatively thick loops to satisfy k > k c . The high density ratio, i.e., no/n e 3> 1 or VAe 3> va, yields the 
following simple expressions for the cutoff wave number k c , 

k c a « j , s (—} = k.s (—\ (5) 



Since the wavelength of the fundamental eigenmode (with harmonic number N = 1 in longitudinal direction 
and mode number s = 1 in transverse direction) corresponds to the double loop length, so that the 
wavenumber relates to the loop length byfc = 27r/A = 7r/Z, the cutoff wave number condition k > k c implies 
a constraint between the loop geometry ratio w/l and the density contrast n e /n n , 

< = \ — ~°- 65 \/— ' ( 6 ) 
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The numerical factor 0.65 applies to the fundamental (s — 1) sausage mode. Since geometric parameters 
(such as I and w) can be measured easier than densities, we might turn the cutoff condition around and 
formulate it as a density contrast requirement for a given loop aspect ratio, 

=■>( > i) 2 = 2 .4(i)\ (7 , 
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This clearly indicates that slender loops with a high length-to- width ratio l/w 1 would require 
extremely high density contrast no/n e . Typical active region loops, which have only a moderate density 
contrast in the order of n /n e w 2— 10, would be required to be extremely fat, with width-to-length ratios of 
l/w s=s 1 — 2. The density contrast is much higher for flare loops or postflare loops, up to no/n e w 10 2 , ...10 3 . 
In this case, a length-to- width ratio of l/w w 6 — 20 would be allowed for sausage mode oscillations. This 
brings us to the conclusion that global sausage type oscillations are only expected in fat and dense loops, 
basically only in flare and postflare loops, a prediction that was not appreciated until recently (Nakariakov 
et al. 2003). The restriction of the sausage mode wave number cutoff is visualized in Fig. 2, where the 
permitted range of geometric loop aspect ratios (l/w) is shown as function of the density contrast (no/n e ). 

Since the sausage mode is highly dispersive in the long- wavelength part of the spectrum (see Fig. 1), 
the phase speed v p h is a strong function of the wavenumber ka. The phase speed diagram Fig. 1 shows that 
the phase speed equals to the external Alfven speed at the long- wavelength cutoff, v p h(k = k c ) = va c , and 
tends to approach the internal Alfven speed in the short- wavelength limit, v p h(ka » 1) i— > va- For coronal 
conditions, the phase speed of both the fast kink and fast sausage mode are bound by the internal and 
external Alfven velocities, 

va < v ph = — < v Ae , (8) 
Therefore, also the period (P = 2l/v p h) of the standing sausage mode is bound by these two limits, 
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At the lower limit we can derive a simple relation for the sausage mode period from the long-wavelength 
cutoff, where v p h(k = k c ) = va £ , using Eq. 6, 
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an approximation that is also derived in Rosenberg (1970) and Roberts ct al. (1984). However, one should 
be aware that this relation represents only a lower limit that applies at k — k c , while for all other valid 
wavenumbers (k > k c ) the loop length is shorter (I — ir/k) and thus the global sausage mode period P aaus 
is also shorter (Nakariakov et al. 2003), given by the actual phase speed v ph (k = ir/l) defined by the 
dispersion relation (Fig. 1) at a particular value k > k c , 

21 

Psaus{l) = v ph (k = * /I) ' (U) 

Provided that sufficiently fat and overdense loops exist according to the requirement (Eq. 6 or 7) , we expect 
for a loops radius of a « 1000 km and an Alfven speed of va ~ 1000 km s _1 a typical sausage mode period 
of Psaus < 2.6 s (according to Eq. 10). 



3. OBSERVATIONS 

Table 1 contains a compilation of radio observations of oscillation events with periods in the range of 
0.5 s < P < 1 min. Most of these observations have been interpreted in terms of MHD fast sausage mode 
oscillations, mainly for three reasons: (1) Eigen-modes of standing MHD waves provide a natural mechanism 
for oscillations with a strictly regular period; (2) the period range of P m 0.5 — 5 s was expected for the 
MHD fast sausage mode under coronal conditions, while other MHD oscillations modes have significantly 
longer periods; and (3) the MHD fast sausage mode modulates the loop density, magnetic field strength, 
and related radio emission (plasma emission or gyrosynchrotron emission), while the MHD kink mode does 
not modulate the loop density in first order. In addition we might add a forth reason which has not been 
appreciated previously: These events listed in Table 1 were generally reported to occur during flares, which 
explains the high densities that are expected on theoretical grounds (due to the sausage mode wave number 
cutoff criterion k c > k, see Eq. 7), which does not occur for the MHD kink mode. 

Because the loop density plays a critical role for MHD fast sausage mode oscillations (Eq. 7), we listed 
in Table 1 only observations that have information on the period P as well on the electron density no of 
the oscillating loops. There are only two observations with direct electron density measurements, using soft 
X-ray emission measures and constraints from the microwave spectrum, which report electron densities of 
n = 4.5 x 10 10 cm" 3 (Asai et al. 2001) and n = 10 11 cm" 3 (Melnikov et al. 2002). However, we might 
also extract information on the electron density from any radio burst that is produced by an emission 
mechanism near the plasma frequency, since the plasma frequency depends only on the electron density n , 

v v = 8980v/n [cm- 3 ] [Hz] . (12) 

Plasma emission is the dominant emission mechanism for metric and decimetric radio bursts up to 
frequencies of v < 1.5 GHz, while free-free emission, gyroresonance, and gyrosynchrotron emission dominate 
at higher frequencies (e.g., see Dulk 1985). A cursory look at Table 1 shows that most of the reported 
oscillations fall in the frequency range of v s=a 100 — 1000 MHz, so this corresponds to electron densities of 
n e ?y 10 8 — 10 10 cm~ 3 in the case of plasma emission. In Fig. 3 we display the observed oscillation periods 
P as function of the observed radio frequency v. 

The observed radio frequency uniquely constrains the density no of the oscillating loop segment, for the 
case of plasma emission. In order to estimate the external density n e we can use a standard model of the 
average background corona, for instance the Baumbach- Allen model, which as function of the normalized 
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height R=(l + h/R Q ) is, 

o T2.99 1.55 0.036 
n e (h) = 10 8 



R ie R e R i.5 



4 x 10 s . 

— ^g- [cm- 3 ] (13) 



where the approximation on the right-hand-side holds for h < Rq/2. This density model n e (h) is shown 
in Fig. 4 (top). Assuming that the fattest possible loops have a width-to- length ratio of q w w 1/4, the 
minimum required density ratio is n ^ min (h)/n e (h) = 2A/q^ w 40 (Eq. 7) for loops oscillating in the MHD 
fast sausage mode due to the long- wavelength cutoff criterion (Eq. 5). This cutoff criterion implies the 
constraint that physical solutions are only possible for heights h £ 40 Mm, for oscillations observed in the 
plasma frequency range of v p < 1 GHz (grey area in Fig. 4 top). The two events observed by Asai et al. 
(2001) and Melnikov et al. (2002) do not fall in this range, because their density is higher and would 
correspond to plasma frequencies of v v = 2 and 3 GHz. The important conclusion that follows, based on 
the two assumptions that the Baumbach- Allen model represents a realistic density model of the background 
corona and that oscillating loops are not fatter than a quarter of their length, is that physical solutions for 
the MHD fast sausage mode are only possible for coronal heights of h £ 40 Mm in the observed frequency 
range of v < 1 GHz, where most of the fast oscillations have been observed (Fig. 3). 

Next we investigate the constraints introduced by the pulse periods, which are found in the range of 
Psaus ~ 0.5 — 5.0 s (Fig. 3). Denoting the length of the loop segment that participates in the MHD fast 
sausage mode oscillation with A/2, we obtain an upper limit from the phase speed at the wavenumber cutoff 
k c (Eqs. 9, 11), 

^ -Psaus (■, \ - Psaus /-* A \ 

2 = ~2~ v ph( k > < — g— V Ae , (14) 
where the Alfvenic speed VAe is given by the background density model n e (h) (Eq. 13), 
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where we used a mean molecular weight of /i = 1.3 for the H:He=10:l coronal abundance. We plot the 
sausage mode period P sa us as function of the wavelength A/2 in Fig. 4 (bottom), for a mean coronal 
magnetic field of B = 20 G (solid line), as well as for a factor of 2 smaller or larger values (B = 10 — 40 
G). We find that the corresponding wavelength of the oscillating loop segments have values in the range 
of 0.3 Mm < A/2 < 10 Mm (gray area in Fig. 4 bottom) for pulse periods of 0.5 < P saU s < 5.0 s, which 
is much smaller than the full loop length (I — irh > 120 Mm, see gray area in Fig. 4 top) required for 
the fundamental sausage mode. Because physical solutions are only possible in heights h > 40 Mm, the 
oscillating loop segment has to be located near the top of the loop. Thus we conclude that only a small 
segment of the loop is oscillating in the MHD sausage mode. Such small oscillating loop segments require 
wavelengths that correspond to higher harmonics N of the fast sausage MHD mode, 

\ 2w 21 fin 

a = T = a' (16) 

while lower harmonics cannot oscillate due to the wave number cutoff. The situation is illustrated 
for the case N = 3 in Fig. 5: For instance, if the cutoff is k > k c — 6n/l, only wavelengths with 
A = 2ir/k < 2ir/k c = 1/3 are able to oscillate, which is most likely to occur at the loop top where the density 
contrast and wavelength match the sausage wavenumber cutoff criterion. Alternatively, oscillations at a 
higher harmonic could occur along the entire loop, but favorable conditions for detection could be restricted 
to certain segments only (e.g., due to line-of-sight angles or absorption of plasma emission). 
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An alternative interpretation of these fast radio pulsations in the decimetric range could be made in 
terms of propagating fast sausage mode MHD waves, where the wavelength is not prescribed by the loop 
length, but independently by the driver. This interpretation has been proposed in Roberts et al. (1984). 
However, if this would be the case, the radio dynamic spectra should show parallel drifting bands of 
emission for every propagating node. Such decimetric burst types have indeed been observed and are called 
decimetric Eber bursts (e.g., Rosenberg 1972; Bernold 1980; Slottje 1981), but many of the observations 
listed in Table 1 show narrowband pulsations without parallel drifting bands, and thus are not consistent 
with (impulsively generated) propagating fast sausage MHD waves. 

Note that the measurements of Asai et al. (2001) and Melnikov ct al. (2002) have longer periods and 
higher densities than all decimetric observations listed in Table 1, and thus actually allow for global sausage 
mode oscillations of the entire loop length, for reasonable magnetic fields in the order of B s=a 40 G (marked 
with crosses in Fig. 4 bottom). 

4. CONCLUSIONS 

This study leads us to the following conclusions: 

1. Only the two imaging observations of fast oscillation events (P « 10 s) by Asai et al. (2001) and 
Melnikov et al. (2002) are consistent with global MHD fast sausage mode oscillations, satisfying 
the relation between period and loop lengths as well as the wave number cutoff criterion for the 
fundamental sausage mode. The required magnetic field is of the order of 5 « 40 G. 

2. All earlier radio observations of fast oscillation events, which have been reported in the frequency 
range of v =100 MHz — 1 GHz and with periods of P ~ 0.5 — 5.0 s, imply for reasonable coronal 
magnetic field strengths (B s=a 10 — 40 G) that the length of the loop segments participating in 
the MHD fast sausage mode oscillation is much smaller (A/2 < 10 Mm) than the full loop length 

(I 100 — 500 Mm). The sausage mode oscillation seems to be confined to segments where the largest 
density contrast to the background occurs, most likely near the loop top, in a segment corresponding 
to a higher harmonic node (see Fig. 5). 

3. The wave number cutoff criterion seems to introduce a bimodality of MHD fast mode sausage 
oscillation phenomena: (1) partial loop oscillations confined to loop segments (A/2 <C I) of high 
harmonic (N ^> 1) nodes for low-density (n e < 10 10 cm~ 3 ) loops that are observable at the plasma 
frequency in the metric/decimetric frequency range [y < 1 GHz); and (2) global loop oscillations 
(A/2 = I) at the fundamental mode (N = 1) for high-density (n e 10 10 — 10 11 cm~ 3 ) loops that are 
observable in soft X-rays and via gyrosynchrotron emission in microwaves. 

In this study we explored the diagnostic of MHD fast sausage mode oscillations. It is important to notice 
that the observed period does not provide a diagnostic of the Alfvcnic travel time across the loop diameter 
(P = 2.62a/ v a), as it was commonly believed, but rather provides a lower limit of this quantity (Nakariakov 
et al. 2003). Moreover, the new finding that fast periods are likely to be caused by oscillations of partial 
loop segments (at higher harmonic nodes), implies that the magnetic field B cannot be determined by the 
loop diameter and sausage mode period, unless the harmonic number can be constrained from imaging 
observations. This study illustrates specific constraints that are contingent on the sausage wave number 
cutoff, which have been ignored in previous studies. 
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Table 1. Observations of oscillations in the period range of 0.5 s < P < 1 min with density information 
(either from radio frequencies of v < 1.5 GHz or from soft X-rays emission measures). 



Observer 


Number 


Frequency 


Period 


Instrument / 




of events 


v [MHz] 


P [s] 


Spatial resolution 


Dr5gc (1967) 


> 18 


240, 460 MHz 


0.2-1.2 


Kiel 


Rosenberg (1970) 


1 


220-320 MHz 


1.0-3.0 


Utrecht 


De Groot (1970) 


> 25 


250-320 MHz 


2.2-3.5 


Utrecht 


Abrami (1970, 1972) 


3 


239 MHz 


1.7-3.1 


Trieste 


McLean ct al. (1971) 


1 


100-200 MHz 


2.5-2.7 


Culgoora 


Rosenberg (1972) 


1 


220-320 MHz 


0.7-0.8 


Utrecht 


Gotwols (1972) 


1 


565-1000 MHz 


0.5 


Silver Spring 


Kai & Takayanagi (1973) 


1 


160 MHz 


<1.0 


17' (Nobeyama) 


McLean & Sheridan (1973) 


1 


200-300 MHz 


4.28±0.01 


Culgoora 


Achong (1974) 


1 


18-28 MHz 


4.0 


Kingston 


Tapping (1978) 


14 


140 MHz 


0.06-5 


Cranlcigh 


Pick &z Trottct (1978) 


1 


169 MHz 


0.37, 1.7 


5'-7' (Nancay) 


Elgaroy (1980) 


8 


310-340 MHz 


1.1 


Oslo 


Bernold (1980) 


> 13 


100-1000 MHz 


0.5-5 


Zurich 


Slottjc (1981) 


> 40 


160-320 MHz 


0.2-5.5 


Dwingcloo 


Trottet et al. (1981) 


1 


140-259 MHz 


1.7±0.5 


5' (Nancay) 


Sastry ct al. (1981) 


1 


25-35 MHz 


2-5 


30' (Gauribidanur) 


Zlobec et al. (1992) 


1 


333 MHz 


9.8-14.2 


0.7'-1.5' (VLA) 


Zaitsev et al. (1984) 


23 


45-230 MHz 


0.3-5 


Izmiran 


Wichl ct al. (1985) 


1 


300-1000 MHz 


1-2 


Zurich 


Aschwanden (1986, 1987b) 


10,60 


300-1100 MHz 


0.4-1.4 


Zurich 


Aschwanden & Benz (1986) 


10 


237, 650 MHz 


0.5-1.5 


Zurich 


Kurths & Herzcl (1986) 


1 


480-800 MHz 


1.0, 3.5 


Trcmsdorf 


Aurass ct al. (1987) 


1 


234 MHz 


0.25-2 


Trcmsdorf 


Kurths & Karlicky (1989) 


1 


234 MHz 


1.3, 1.5 


Tremsdorf 


Chernov & Kurths (1990) 


10 


224-245 MHz 


0.35-1.3 


Izmiran 


Kurths ct al. (1991) 


25 


234-914 MHz 


0.07-5.0 


Zurich 


Aschwanden ct al. (1992c) 


1 


1.5 GHz 


8.8 


0.2'-0.9' (VLA, OVRO) 


Aschwanden et al. (1994a. b) 


1,1 


300-650 MHz 


1.15, 1.8 


Zurich 


Chernov ct al. (1998) 


1 


164-407 MHz 


0.2 


5'-7'(Nancay) 


Kliem ct al. (2000) 


1 


600-2000 MHz 


0.5-3.0 


Zurich 


Asai ct al. (2001) 


1 


17 GHz 


6.6 


10" (Nobcyama) a 


Melnikov ct al. (2002) 


1 


17, 34 GHz 


8-11, 14-17 


10", 5" (Nobcyama) fa 


a Asai et al. (2001) estimated 


a loop length of L = 16 Mm, a loop 


width of in = 6 Mm, 


and a loop density of n e = 4.5 X 


10 10 cm -3 , using Nobeyama and 


Yohkoh/SXT data observed during the 1998-Nov-10, 00:12 UT, flare. 








b Melnikov ct al. (2002) estiir. 


^ated a loop length of L = 25 Mm, a 


loop width of w = 6 


Mm, and a loop density of n e = 


10 11 cm -3 , using JVobcyama and 


Yohkoh/SXT observed during the 2000-Jan-12, 01:35 UT. flare. 
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Fig. 1. — The phase speed u/k is shown for magnetoacoustic waves in a cylindrical fhrxtube (with radius 
a), as function of the longitudinal wave number ka, for coronal conditions VAe > va > Ct > c s . The sausage 
modes are indicated with solid lines, the kink modes with dashed lines (Edwin & Roberts 1983). 
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Fig. 2. — The wave number cutoff k c for sausage mode oscillations expressed as a requirement of the loop 
length-to- width ratio I /w as function of the density contrast n /n e between the external and internal loop 
densities (see Eq. 6). Note that sausage mode oscillations in slender loops can occur only for flare conditions. 
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Fig. 3. — Observed oscillation periods as function of the radio frequency according to the list given in Table 
1. Horizontal and vertical bars indicate the ranges in each observation, diamonds the mean values, and 
thick diamonds label reports on multiple oscillation events. Note that most of the events fall in the range 
of P m 0.5 — 5.0 s and n = 10 8 — 10 10 cm~ 3 , if one assumes that the radio frequency is near the plasma 
frequency v w v p . The letter A and M refer to the observations of Asai et al. (2001) and Mclnikov ct al. 
(2003). 
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Fig. 4. — Top: The Baumbach-Allen density is shown as a model of the average coronal background (thin 
curve), together with the approximation n e (h) «4x 10 8 (1 + h/R Q ) 9 (thick curve), and the minimum loop 
density n ^ min (h) required for loops that have a width-to-length ratio of q w = w/l < 0.25. The height 
and density range corresponding to a plasma frequency range of 100 MHz-1 GHz is indicated with a grey 
area. Bottom: The MHD fast sausage mode period P saU s is shown as function of the sausage length A/2 for 
magnetic fields of B = 20 G within a factor of 2. The regime of physical solutions for periods P — 0.5 — 5 s 
and B = 10 — 40 G is indicated with a grey area. A and M indicate the values inferred for the observations 
of Asai et al. (2001) and Melnikov et al. (2002). 



- 13 - 




Fig. 5. Sketch of MHD fast sausage mode oscillation in partial loop segment (right), where the cutoff 
condition k > k c is satisfied. For a given width-to-length ratio q w , oscillations occur in segments where 
no/n e > (Eq. 7). The oscillating segment may correspond to a higher harmonic node, e.g., to N = 3 

here (left). 



